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Torsion models constitute a well known class of extended quantum gravity models. In this paper 
we study some phenomenological consequences of a torsion field interacting with fermions at LHC. A 
torsion field could appear as a new heavy state characterized by its mass and couplings to fermions. 
These new states will form a resonance decaying into difermions, as occurs in many extensions of 
the Standard Model, such as models predicting the existence of additional neutral gauge bosons, 
usually named Z' . Using the dielectron channel we evaluate the integrated luminosity needed for 
a 5a discovery as a function of the torsion mass, for different coupling values. We also calculate 
the luminosity needed for discriminate, with 95% C.L., the two possible different torsion natures. 
Finally, we show that the observed signal coming from the torsion field could be distinguished from 
a signal coming from a new neutral gauge boson, provided there is enough luminosity. 

PACS numbers: 13.85.-t,04.60.-m 



I. INTRODUCTION 



One great challenge of the Standard Model is the diffi- 
culty to incorporate gravity. Several alternative theories 
of quantum gravity, extending general relativity, intro- 
duce an extra field called torsion [H, H, Ej] , where the spin 
of the elementary particles is the source of the torsion. In 
a phenomenological point of view, torsion can be treated 
as a fundamental propagating field, either of vectorial or 
pscudo-vectorial nature, characterized by its mass and 
the coupling values (77/) between torsion and fermions. 
Despite there are strong theoretical restrictions on the 
possible parameters of the propagating torsion field 0], 
there is still room for detecting this field on high en- 
ergy colliders. The introduction of torsion fields leads to 
singularity-free cosmological models 0, @] and they can 
be investigated in the near future high energy particle 
colliders, such as the LHC at CERN, through their inter- 
action with the fermions. Some torsion studies have been 
previously performed for LEP and TEVATRON @, @] 
and even for the LHC Other torsion aspects have 
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been also recently studied 13, U|- The LHC, which will 
have a center-of-mass energy of 14 TeV, will provide un- 
precedented conditions to investigate many alternative 
models and will extend the search for new heavy parti- 
cles up to the 5-6 TeV range, depending on the integrated 
luminosity obtained. The machine must start its activity 
at low luminosity (C = 10~ 33 cm^s^ 1 ), reaching high 
luminosity (C — 1CP 34 cm -2 s -1 ) after one year of opera- 
tion 

Here we present the discovery potential of torsion field 
at the LHC. The dielectron channel is used to calculate 
the luminosity needed for a 5er discovery as a function of 
the torsion mass and the fermion-torsion vertex coupling. 
We consider all the couplings 77/ equal. The analysis with 
early LHC data of new heavy states formed by opposite 
sign dileptons will have great importance, since these sig- 
natures are included in many Standard Model extensions. 
The LHC experimental data can eliminate or corroborate 
some of these proposed phenomenological extensions. 

In a second part of this work, we present the luminosity 
needed for the discrimination of the torsion field nature, 
vectorial or pseudo- vectorial, with 95% C.L.. Finally, we 
present a variable capable of separating the torsion field 
from a new vector neutral gauge boson (Z'), predicted 
by several models, such as the E§ [ijj]. We show that, 
once a signal of a dilepton resonance is experimentally 
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found, the accumulated asymmetry (defined in Sec. V) 
can discriminate between the torsion and the Z' . The 
parameters and particle efficiencies of the ATLAS detec- 
tor have been roughly used in this paper, although no 
detector simulation had been performed. 

II. THE TORSION FIELD 

There are many experimental and theoretical argu- 
ments saying that the minimal Standard Model (SM) is 
not a complete theory and it should be extended. One 
possibility is to extend it in such way to include also grav- 
ity. This inclusion should predict the observable value of 
the cosmological constant and bring new observable ef- 
fects in high energy elementary particle physics. The 
SM is described using three types of fields: scalar, vec- 
torial and spinorial, by the other hand general relativity 
has in addition the metric field which describes the geo- 
metric property of space-time. A strong candidate to be 
included in the SM is the space-time torsion field which 
adds some independent characteristics to the space-time 
geometry as is shown in the general relativity with tor- 
sion [2]. The torsion field can be considered as being 
composed of three irreducible components: the axial vec- 
tor S^, the vector trace T a , and the tensor 

The general minimal and non minimal action of a Dirac 
fcrmion coupled to torsion can be written for practical 
purposes as 



(1) 



S f = I d*XT/g{ijn»(y ll -ir) 1 >y* i S li 



where r\\ and 772 are two dimensionlcss parameters and 
V« is the Riemannian covariant derivative (without tor- 
sion). The q a p^ tensor decouples completely from the 
fcrmion fields for the minimal and non minimal actions. 

One has two approaches: the minimal and the non 
minimal interaction. The minimal interaction gives 771 = 
— 1/8 and the vectorial torsion part decouples completely 
from the spinor fields. This case is equivalent to put 
772 = in Eq. 1 and then one has a pure axial vector 
coupling. So for the minimal case the theory gives a 
fixed value of 771 which corresponds to a pure axial vector 
coupling [2j . 

By the other hand, the non minimal interaction allows 
both rji and 772 to have non zero values. In order to 
simplify our study and analyze the extreme cases, we are 



considering the pure axial vector (77/ = 771 7^ and 772 =0) 
and the pure vectorial (77! =0 and r\f — 772 7^ 0) cases, 
although the mixed case (771 7^ and 772 7^ 0) is not at 
all excluded. We analyze the non minimal case since the 
minimal one is very restrictive due to the fixed values of 
771 and 772 in Eq. 1. 

For the axial case, the non minimal interaction of 
torsion with Standard Model fermions will be given by 
the following action, where ip^ stands for each of the 
fermions: 



^TS -matter _ ■ I A 



+imH 5 7 tl S ll - irrii ^>(i) (2) 

While for the vectorial case, the non minimal action 
will be: 



pT 'S— matter 
V—non—min 



^Vg (7" v Q 



(3) 



where 77H and 772^ are the interaction parameters for the 
corresponding fcrmion spinor and the subscripts V and A 
stand for the vectorial and pseudo-vectorial cases. From 
the theoretical point of view the nonminimal interaction 
with the geometric vector component of torsion T M can 
not be ruled out [2] and therefore it is worthwhile to 
work out it in the same manner as it is done with the 
axial vector component. 

We have considered the non minimal parameters equal 
for all fermions, contrary of what was done on [9( where 
the top quark coupling was taken to be different from the 
others due to its larger mass because of the renormaliza- 
tion group running. So, from now on, 777- = r\u = if 
r/2 = and 77/ = 772; = 772, if 771 =0 then we are dealing 
with two parameters only, 77/ and Mts for the axial and 
vectorial cases. We are studying the channel pp — > e + e~ 
and this assumption with respect to the top quark cou- 
pling makes practically no difference in the results here 
obtained. This assumption changes slightly the torsion 
width. 



III. EVENT SELECTION 

Torsion (TS) described as pure vectorial and pure 
pseudo-vectorial field was implemented in CompHep 
event generator [l4|. Using this package, we simulate the 
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process pp — > j/Z/TS — > e + e _ with a center-of-mass 
energy of 14 TeV and with full interference between the 
torsion and the SM gauge bosons. The events were gen- 
erated for three different values of the torsion mass: 1, 
3, and 5 TeV and for three values of rjf. 0.2, 0.5 and 0.8, 
which obeys the limits obtained by 0] and the constraint 
Mrs/Vf ^> m/ (or, equivalently, Mts ^ Vf x m /)? ob- 
tained in (jj for the pure axial case. For the pure T M case, 
there are no such theoretical limits. A cut on the dielec- 
tron invariant mass of 500 GeV was applied on gener- 
ated events in order to reduce the SM background. The 
CTEQ6L parton distribution functions were employed. 
The pseudo-rapidities of the electron and the positron 
were required to be \rj\ < 2.5 and each final particle en- 
ergy should be greater than 5 GeV, both conditions due 
to the ATLAS detector characteristics. Since in a more 
realistic scenario jets can be misidentified as electrons in 
calorimeters, which affect the electron identification, we 
assume a dielectron identification efficiency of 62%. The 
overall efficiency, which take into account the cut on r\ 
and the electron identification is around 50% New 
neutral gauge models based on Eq were also implemented 
in CompHep. 



IV. TORSION DISCOVERY POTENTIAL AT 
LHC 

In Fig. [TJ we present the torsion width as a function of 
its mass for several fermion-torsion couplings. Due to its 
huge width, in special for large rjf values, torsion inter- 
action can not be approximated as a contact interaction. 
Figure [2] shows the production cross section at NNLO for 
the process pp — > j/Z/TS — > e + e~ as a function of tor- 
sion mass and for three different values of rjf . LO cross 
section obtained with the CompHep package was multi- 
plied by a constant factor K = 1.35 in order to take into 
account the next-to-next-leading order QCD corrections. 
As expected, the cross section grows as ?7/ increases, since 
the fermion-torsion coupling becomes stronger. We can 
see that for Mts above 4 TeV, there is no significant dif- 
ference between the cross sections for the different cou- 
plings, showing clearly the effect of the parton distribu- 
tion functions. The widths and total cross-sections are 
equal for both the pure axial and the pure vectorial cases. 

The dielectron channel was chosen to calculate the 
torsion discovery potential at LHC since ATLAS has a 
better electron resolution when compared to muon 
A signal of a torsion field interacting with the fermions 
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FIG. 1: Torsion widths as a function of torsion mass for sev- 
eral fermion-torsion couplings (rjf = 0.2,0.5 and 0.8). 
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FIG. 2: Torsion production cross section as a function of tor- 
sion mass, for several fermion-torsion couplings (77/ = 0.2, 0.5 
and 0.8). 



should be observed above the neutral Drell-Yan (DY) 
process that is the main background in the search of new 
heavy dilepton resonances. The statistical significance of 
the torsion signal in the presence of background can be 
obtained using the likelihood-ratio estimator 



S L = ^2\n{L s+b /L b ) 



(4) 



where L s+b represents the maximum likelihood value ob- 
tained from the unbinned likelihood fit to the dielectron 
invariant mass assuming that a signal is presented in the 
data, while L b is the maximum likelihood value obtained 
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FIG. 3: Background shape determination for M ee < 3 TeV 
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FIG. 4: Background shape determination for M ee > 3 TeV 



with a background only hypothesis. It has been shown 
that this estimator has the desired feature of not having 
an overestimated or underestimated probability of a false 
discovery [l7| . 

The probability density function (pdf) used to fit 
the observed dielectron invariant mass spectra is given 
by Q. 



TABLE I: Cuts on M ee obtained from the exponential fit to 
the invariant mass for different values of Mrs and rjf. 





Torsion Mass (Mrs) 


Vf 


ITeV 


2 TeV 


5 TeV 


0.2 


660 


1570 


3500 


0.5 


680 


1360 


3300 


0.8 


680 


1400 


3300 



p (M ee ) = f sPs + (1 - fs)p b (M ee ) (5) 

where p s is the pdf of the signal, given by a Breit-Wigner 
distribution, and pb is the pdf of the background, mod- 
eled as exponential functions. The parameters to be fit- 
ted are the signal fractions f s — N S /(N S + Nb) (where 
Ngn,) are the number of signal (background) events), the 
torsion mass Mts peak position and the torsion width 

The pdfs that describe the background shape depends 
on the invariant mass region, and are given by 

(expi-hM^ 2 ), if 500 < M ee < 3 TeV 
(6) 
exp(-k 2 M ee ), if M ee > 3 TeV 

The values of the parameters k\ and ho. are obtained 
from fits to the DY invariant mass Monte Carlo (MC) 
distribution in the full mass region of interest. Figures 
12] and 2] show the fits and the resultant values of k\ and 
k 2 : h = 6.38 and k 2 = 0.00209. We can see that the 
DY background can be well described by the proposed 
parametrization. 



In order to increase the signal-to-background ratio, it 
was applied a M ee cut to eliminate as much as possible 
the background contamination on the mass signal. Fig- 
ures [5] and [6] illustrate for two different values of torsion 
masses how the cut is obtained: the Equation 6 is used to 
fit the invariant mass distribution, with fex,2 fixed at the 
values mentioned above, and with the amplitude been 
the only fitted parameter. The M ee cut is defined as the 
point where the da/dM ee reaches the value 0.5/6. Table 
U shows the cut values obtained for the different torsion 
masses and couplings considered in this study. 

To obtain the statistical significance (Eq.4), Equation 
5 is fitted to the invariant mass of the events that pass the 
cut described above. The fit is performed in two steps: 
first, using Equation 5, we fit the invariant mass for many 
different MC samples, and take the mean values of Mts 
and Tts from the fitted parameters distribution. Then, 
fixing Mts and Tts at their mean values, we fit the 
same samples again, under the signal plus background 
hypothesis (note that in this case, just f s is fitted), in 
order to obtain L s+ b- By making f s = in Equation 
5, which represents the background only hypothesis, we 
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FIG. 5: Mass cut determination for Mrs = 3 TeV 
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FIG. 6: Mass cut determination for Mrs = 5 TeV 
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FIG. 7: Significance distribution Sl for Mrs = 1 TeV and 
W = 0.2 
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fit the same distributions to get In the last fit the 
parameters k± and &2 are free. As the result of many fits 
a Sl distribution is obtained, as illustrated for the case 
Mrs = 1 TeV and rjf = 0.2 in Figure [7) The mean value 
of Sl distribution is taken as the signal significance for a 
given luminosity. The value of Sl must be greater than 
5 for claiming a discovery. 

Fig. [5] shows the integrated luminosity needed for a ha 
discovery as a function of torsion mass for three different 
values of torsion-fermion couplings. From this plot, we 
obtain the following conclusions: i) a very low integrated 
luminosity, of order of 100pb _1 , is needed for discovering 
a torsion of 1 TeV mass. It means that such signal could 
be observed at very early days of LHC operation, even 
in a regime of low luminosity; ii) for a torsion mass of 3 



FIG. 8: Minimum integrated luminosity needed for a 5 a 
discovery 



TeV, the discovery would be possible with an integrated 
luminosity around 10 fb -1 , or approximately 1 year of 
LHC operation in low luminosity; iii) finally, around 10 
years of operation in high luminosity would be needed 
for the observation of a 5 TeV mass resonance. How- 
ever, for the case Mts — 5 TeV and r\j = 0.8, no peak 
can be observed at LHC due to the huge torsion width 
(Tts = 1-9 TeV), although some excess above the DY 
background can still be identified. For this reason, this 
particular point is not included in Fig. [H] 
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V. TORSION NATURE DETERMINATION 

We obtain now the minimum integrated luminosity 
needed to discriminate between the axial and the vec- 
torial models with 95% C. L. as a function of the torsion 
mass for different ry/ values. Assuming that a torsion 
signal was found at LHC and its mass and width were 
already experimentally determined, we have considered 
a mass window for the dielectron invariant mass between 
Mts/2 and Mts — IVs, where Mts and Tts are the 
fitted values for the torsion mass and the torsion width, 
respectively. This is the region where the interference 
terms dominate the process. 

The pseudorapidity difference An between the electron 
and the positron can be used to distinguish between the 
axial and the vectorial torsion models. Due to the en- 
ergy cuts on the final electron energies (i? e > 5GeV) An 
is Lorentz invariant. Figj9] shows the An distribution 
for Mts — 3 TeV and rjf — 0.2 for an integrated lu- 
minosity of 165 fb^ 1 . In order to quantify and test the 
statistical compatibility of the two histograms and since 
either can have bins with few or even zero contents, we 
use the % 2 — function fl!^ of An distributions to cal- 
culate the minimum integrated luminosity that yields a 
95% C.L. discrimination between the axial and vectorial 
models. FigflOl shows our results as a function of the 
torsion mass for 77/ = 0.2, 0.5 and 0.8. We can see that 
for Mts — 1 TeV and r/f = 0.8, an integrated luminos- 
ity of 0.2 fb _1 , corresponding to less than one day of 
LHC data-taking at high luminosity, will be enough to 
discriminate between axial and vector torsions, while for 
Mts = 5 TeV and rjf = 0.2, an integrated luminosity of 
1500 fb _1 , corresponding to more than 16 years of LHC 
data-taking at high luminosity, will be needed to achieve 
the torsion nature distinction. 



VI. TORSION AND NEW NEUTRAL GAUGE 
BOSON DISCRIMINATION 




FIG. 9: Pseudorapidity difference between the electron and 
the positron for axial (solid line) and vectorial (dotted) mod- 
els, for Mts = 3 TeV and nj = 0.2 for an integrated lumi- 
nosity of 165 fb~ . 




FIG. 10: Minimum integrated luminosity needed for 95% C.L. 
discrimination between the axial and the vectorial models us- 
ing the electron-positron pseudorapidity difference distribu- 
tion. 



As the torsion field is characterized by a neutral gauge 
boson, its mass signal at LHC could be the same of a new 
heavy neutral resonances, like Z' . In order to identify the 
signal nature, if its is found, we define an accumulated 
asymmetry, that is given by the forward-backward asym- 
metry as a function of M ee evaluated at increasing ranges 
of the final electron invariant mass M ee . The lower range 
limit is at a fixed value M-rs/2 and gradually increasing 
the upper limit up to Mts + 500 GeV. 



A(M ee ) = 



Np — N 



B 



N F + N B 



M, 



(7) 



where Np is the number of forward events, counted as the 
number of events with cosO* > 0, and Nb is the number 
of backward events, which have cosO* < 0, where 9* is 
the angle between the outcoming electron and the boost 
direction at the center of mass frame of the final elec- 
trons. The errors of the accumulated asymmetry become 
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FIG. 11: Discrimination between the axial torsion, vectorial 
torsion and Z using the accumulated asymmetry defined in 
Equation 7 for Mrs = M z , = ITeV and rjf = 0.2. 



smaller as the invariant mass range increases since the 
total numbers of events also increases. In order to com- 
pare the Z' and torsion fields, we have considered the 
Z' model, but the analysis can be equally applied to any 
neutral gauge boson coming from different models. 

Figure [TT] shows the accumulated asymmetry for the 
torsion fields (vectorial and pseudo- vectorial) and for Z' , 
calculated for the case Mts = Mz' = ITeV, and n/ = 
0.2 with an integrated luminosity of 20fb _1 . We can see 
that the Z' can be easily discriminated from the pseudo- 
vectorial torsion in all range of mass considered. 

On the other hand, the Z' curve overlaps the vectorial 
torsion up to an invariant mass of 1 TeV, but can also be 
distinguished after that point. Note that it is also possi- 
ble to distinguish between the vector and pseudo vector 
torsion fields using this variable. The discriminations are 
easier for greater values of rjf, since the torsion produc- 
tion cross section increases and therefore the errors are 
smaller. 

In Figures [T2] and [T3] we see the same distribution for 
gauge bosons mass of 3 TeV and for rjf — 0.2 and 0.5, 
respectively. In this luminosity of 900 fb is 

needed for discrimination between the resonances. For 
rjf = 0.2, the accumulated asymmetry can distinguish Z 1 
and pseudo-vectorial torsion, but has no power to sepa- 
rate the Z' and vectorial torsion. However, this variable 
shows a good performance on discrimination if 77^ = 0.5, 
as exhibited in Figure [T51 



FIG. 12: Same as Fig. 11 but for Mts = M Z ' = 3 TeV. 
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FIG. 13: Same as Fig. 11 but for Mrs = M z < = 3 TeV and 
V/ = 0.5 



VII. CONCLUSIONS 

We have studied the phenomenological consequences of 
a pure axial vector or pure vector torsion action at CERN 
LHC. In this context, torsion is treated as a massive spin 
1 particle, that will form a resonance and will decay into 
the SM fermions and the model parameters will be the 
torsion mass (Mts) an d the non minimal torsion-fermion 
couplings (j]f). Using the dielectron channel we have cal- 
culated the torsion 5cr discovery at LHC and showed that 
if such a resonance exists, LHC will be able to detect it 
after data-taking corresponding to 100 plV 1 , if the tor- 
sion mass is around 1 TeV or lOflV 1 if the torsion mass 
is around 3 TeV. Unfortunately, for a torsion mass of 5 
TeV and rjf = 0.8, no peak should be observed due to 
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torsion huge width. 

The initial main focus of LHC will be on discovery, but 
after this initial phase, focus will turn to model discrimi- 
nation. Concerning the torsion model, the discrimination 
will have two main points: the determination of torsion 
nature (axial vector or vectorial) and the distinction of 
the torsion field from other new particles predicted by 
many SM extensions, such the Z' , that will also appear 
in proton-proton collisions as resonances decaying into 
difermions. We showed that the pseudorapidity differ- 
ence of the final electron and positron can be used to 
distinguish between an axial and a vectorial torsion. As 
an example, the A77 will be able to discriminate with 95% 
C.L. between axial and vector torsions for Mrs = 1 TeV 
and 77/ = 0.8, with an integrated luminosity of 0.2 fb _1 . 



Another variable capable of doing this distinction is the 
accumulated asymmetry, given by the forward-backward 
asymmetry value evaluated at different mass ranges. But 
this asymmetry will require a greater integrated luminos- 
ity to achieve the distinction between different torsion 
natures and is better suited in the discrimination of the 
torsion fields from an additional neutral gauge boson pre- 
dicted by the Z' x model. When the experimental data be 
available, combining the two methods will certainly in- 
crease the discrimination power of such proposed models. 
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